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(Dated:) 

We consider the ferromagnetic Ising model on a highly inhomogeneous network created by 
a growth process. We find that the phase transition in this system is characterised by the 
Berezinskii-Kosterlitz-Thouless singularity, although critical fluctuations are absent and the mean- 
field description is exact. Below this infinite order transition, the magnetization behaves as 
exp(— const/ \/T c — T). We show that the critical point separates the phase with the power-law 
distribution of the linear response to a local field and the phase where this distribution rapidly 
decreases. We suggest that this phase transition occurs in a wide range of cooperative models with 
a strong infinite-range inhomogeneity. 

Note added. — After this paper had been published, we have learnt that the infinite order phase transi- 
tion in the effective model we arrived at was discovered by O. Costin, R.D. Costin and CP. Griinfeld 
in 1990. This phase transition was considered in the following papers: 

[1] O. Costin, R.D. Costin and CP. Griinfeld, Infinite-order phase transition in a classical spin 
system, J. Stat. Phys. 59, 1531 (1990); 

[2] O. Costin and R.D. Costin, Limit probability distributions for an infinite-order phase transition 
model, J. Stat. Phys. 64, 193 (1991); 

[3] M. Bundaru and CP. Griinfeld, On a phase transition in a one-dimensional non-homogeneous 
model, J. Phys. A 32, 875 (1999); 

[4] S. Romano, Computer simulation study of one-dimensional lattice spin models with long-range 
inhomogeneous interactions, Mod. Phys. Lett. B 9, 1447 (1995). 

We would like to note that Costin, Costin and Griinfeld treated this model as a one-dimensional 
inhomogeneous system. We have arrived at the same model as a one-replica ansatz for a random 
growing network where expected to find a phase transition of this sort based on earlier results for 
random networks (see the text). We have also obtained the distribution of the linear response to a 
local field, which characterises correlations in this system. We thank O. Costin and S. Romano for 
indicating these publications of 90s. 
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The ferromagnetic Ising model on lattices has an ordi- 
nary second-order phase transition [5j- Above the upper 
critical dimension of the model, the critical fluctuations 
are absent, and the mean-field description of this transi- 
tion is exact. In particular, this takes place if couplings 
between all spins are equal — infinite-range interactions. 
In this Letter we report the finding of a phase transition 
with the Berezinskii-Kosterlitz-Thouless (BKT) singu- 
larity in the Ising model on a growing network, which is 
infinitely-dimensional as most of networks. This transi- 
tion is quite unusual for an infinitely-dimensional system 
as well for a cooperative model with the order parameter 
of discrete symmetry. 

Recall that in "ordinary" continuous phase transi- 
tions, pair correlations of an order parameter show a 
slow, power-law space decay only at the critical point, 
T c , and decay exponentially both in the low- and high- 
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temperature phases. This behavior was observed in the 
Ising model on equilibrium complex networks 0, 

Hill 

10] (for percolation and for disease spreading on equilib- 
rium networks, see Refs. ^jj and 113 . respectively). In 
contrast, the BKT phase transition |l3tll4| separates the 
phase with rapidly decreasing correlations and the criti- 
cal phase with correlations decaying by a power law. The 
contact of these two phases is characterised by specific de- 
pendences. For example, the order parameter behaves as 
M(T) ~ exp(— const />/T c — T), and the phase transition 
is of infinite order. 

Normally, the BKT transition is realized in systems 
with two-component order parameters of continuous 
symmetry at a lower critical dimension. Also, this 
anomalous phase transition is present in a few low- 
dimensional systems (e.g., the Luttinger liquid) which 
actually can be reduced to above indicated ones. There 
is one more interesting situation, where the BKT sin- 
gularity emerges. It was observed that in some grow- 
ing networks, near the birth point of the giant connected 
component, its relative size behaves similarly to the mag- 
netization near the BKT transition pi 111 Il7l llH llfll El 
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The model. — We find an exact solution of the following 
cooperative model. A network grows up to a large size, 
and and interacting spins are considered on the resulting 
net. The interaction between the nearest neighbor spins 
is described by the ferromagnetic Ising model. 

We use the following growing network: 

(i) The growth starts with a single vertex (t = 0). 

(ii) At each time step, we add a new vertex and attach 
it to one of "older" vertices. 

(iii) For simplicity, we suggest a specific annealing. For 
an edge born at time t, the end of the edge at vertex 
t is fixed, and the second end can be found at each 
of the vertices in the range < r < t with equal 
probability. Characteristic times for the jumps of 
this end between the vertices < r < t are assumed 
to be not greater than those of the spin relaxation. 

One can show that the resulting model is equivalent to 

the ferromagnetic Ising model on the deterministic graph 
shown in Fig. ^ In this system, the spin on a vertex, 
which was born at time i, has equal coupling l/t to each 
of spins on the older vertices. The Hamiltonian of the 
model is: 



0<i<j<t 



J 



(1) 



where spins = ±1, and Hi > is an inhomogeneous 
magnetic field. Actually, we reduce our problem to a sys- 
tem with a strong deterministic infinite-range disorder. 

Mean-field treatment. — Let us first use a mean-field 
ansatz. We will show afterwards that the mean-field so- 
lution is exact. For the sake of brevity, we use the fol- 
lowing simple mean-field treatment. We assume small 
fluctuations of spins from their mean-field values m,: 
SiSj — > niinij + m,*(sj — rrij) + mj(s,i — mf). Substituting 
this relation into Eq. Q gives a linear effective mean- 
field Hamiltonian. With this Hamiltonian, it is easy to 



obtain the partition function Z 
(f3 = 1 /T) and the free energy F 
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FIG. 1: Deterministic system of interacting spins which is 
equivalent to our growing network. The numbers on the edges 
show the values of the Ising couplings between the spins at 
the corresponding vertices. 



Here we assumed that t is large and passed to the con- 
tinuum limit: mi — m(x = i/t). Expression (J2J together 
with the relation m(x) — —(l/t) SF/5H(x) allows us to 
obtain the equation for the mean local magnetization: 



m(x) = tanh <f3 



- [dym(y)+ f — m(y) + H 

x Jo Jx y 



(x) 

1 w 

The exact derivation of the free energy. — The free en- 
ergy can be found exactly. We compare the free energies 
of the network at times t — 1 and t. For brevity, here we 
consider only the homogeneous magnetic field. The form 
of the Hamiltonian results in the relation 



,-PF t (H) _ 



s t =±l 



e -l3F t ^(H+s t /t) e l3Hs t 



(4) 



that is, 



-l3lF t {H)-F t - 1 (H)} 



-f3F t (H)/t 



Yl exp 

s=±l 



,9F t (H) 



dH t 



(3Hs 



(5) 



Here we took into account the fact that at large t, the 
ratio F t (H)/t approaches a t-independent limit. Using 
the relation for the (relative) full magnetization M(H) = 
J dxm(x.H) = —(l/t)dF t (H)/dH, we get the exact 
form of the free energy: 



F = -t P^ 1 ln{2 cosh[/3(_ff + M(H))]} 



(6) 



at t — ► oo. One can check that free energy ex- 
pressions iPJ — the mean-field one — and © — the exact 
expression — coincide. Indeed, substituting Eq. © into 
the relation (J2J and making partial integration, we ar- 
rive at the free energy exactly in form ©. In this sense, 
the mean-field treatment of this problem is exact. 

Analysis of the equation for the magnetization. — Let 
us consider Eq. From this equation, one can see that 
the assumption m(x) ^ at some x immediately leads 
to the following behavior of m(x) near x = 0: 



m(x ~ 0) = 1 — A x 



20 



(7) 



where A depends on [3 and H . If H = 0, this behav- 
ior is realized only in the low-temperature phase, and 
m(x) = above T c . If H > 0, m(x = 0) = 1 at 
any temperature. (We will see that the critical point, 
T c , exists.) On the boundary, Eq. © readily gives 
m(l) =tanh{/3[M + ff(l)]}. 

For finding this profile, it is convenient to pass to a 
differential equation. For brevity, here we assume that 
H = 0. We introduce a new variable, n(z): 



n(— In a;) = 



- f dym(y) + f — m(y) 

1 Jo Jx y 



so 



m(x) = tanhn(— \nx) . 



(8) 



(9) 
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Differentiating Eq. (|SJ and using Eq. 10 gives the second 
order differential equation 



dn{z) d 2 n(z) 
dz dz 2 



/3tanhn(z) 



(10) 



with the boundary conditions: (i) (dn/dz)(z~0) = n(z — 
0) [note that n(z = 0) = (3M = j3 J Q dym(y)} and (ii) 
n(z — > oo) = (3z + const, z is related to x = i/t: z = 
— In £, so < z < oo, where z = corresponds to x = 1. 
Boundary conditions (i) and (ii) follow from definition 
JSJ) and relation J7J, respectively. At each value of j3, 
there is a single solution of Eq. (|f with these boundary 
conditions, which allows one to get M . 

Equation (|f 0f> can be transformed into a first order 
differential equation. For this, we pass from variables 
{t, n(t)} to {n,w(n)}, where w = f3~ 1 (dn/dt). [n varies 
from to co, while w(n) takes values between and 1.] 
This gives the equation 



w——=/3 1 (w — tanh n) 
dn 



(11) 



for w(n) with the following boundary conditions: (i) 
w[n(z = 0)] = p- 1 n(z = 0) [recall that p- 1 n(z = 0) = M] 
and (ii) w(n — > oo) = 1. Here, boundary condition (i) on 
the line w — f3~ 1 n corresponds to that at x = i/t = 1 . 
Asymptotic boundary condition (ii) corresponds to the 
limit i/t — > 0. Knowing w(n) one can easily get m(x). 

The analysis of Eq. (|1H is similar to that of an equa- 
tion of this type in Ref. 16j. At small n, one can substi- 
tute tanh n by n on the right-hand side of Eq. (fTT)) . so we 
have wdw/dn = (3~ l (w — n). The solutions of this equa- 
tion can be presented in an analytical form. A physically 
reasonable non-zero solution must cross the ordinate axis 
(and the w = line) at non- negative w. This solu- 

tion exists if (3 > 1/4. There is a critical point, f3 c = 1/4, 
where the solution is 



,(»,/?= 1/4) = 2n[l -/(»)] 



(12) 



with f(n) satisfying /(n — > 0) — » and the relation: 
ln[n/(n)] + l//(n) = Inc. Here the constant c = 1.554. . . 
ensures that w c (n) <|12fl fits the corresponding solution of 
Eq. Ijlll) which approaches 1 as n — > oo. The form of the 
critical solution at small n indicates the presence of the 
BKT singularity. 

Near T c , the solution of Eq. (|11|) is close to the critical 
one. In this range, the asymptotics of the solution at 
small n satisfies the relation: 



1 \2f3w(n)/n] - 1 

: arctan- v ;/ J 



V4/T^T 



In \/ v? — w{n)n + f3w 2 {n) — const 



(13) 



(/3 > 1/4). This asymptotics and the solution at large n 
can be sewed together (see details in the full version of the 
present work). For obtaining the dependence of the full 
(relative) magnetization on (3 near the critical (3 C = 1/4, 
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FIG. 2: Main thermodynamic quantities versus (3 = 1/T near 
the phase transition (f3 c — 1/4): (a) the (relative) magnetiza- 
tion, (b) the specific heat, and (c) the magnetic susceptibility. 



we use the following procedure, (i) We substitute the 
boundary condition w(n — (3M) — M into relation 113|l . 
After expansion of the arctangent, we obtain the left- 
hand side of the relation below: 



tt/2 



vr/2 



1 - In 



M{(3) 



w(n) 



2 n 



hi 



w(n) 



tt/2 



In, 



(ii) On the other hand, near f3 ~ 1 /4, in the region (3M <C 
77 <C 1, the main contribution of Eq. I|13|) gives the right- 
hand side of the equality above. We also use the fact that 
the solution must approach the critical one as (3 — » 1/4. 
So, we obtain the full magnetization near (3 C — 1/4: 



M((3) = 4ceexp - 



2^/3 - l/4 y 



(14) 



where 4ce = 16.90 . . ., e is Euler's number. Note that this 
BKT behavior is a direct result of the specific singular 
form of Eq. (|llfl at small n and w. The behaviors of the 
magnetization and other main thermodynamic quantities 
near the phase transition are shown in Fig. [21 

By using Eq. I)ll|l . one can also find the coefficient of 
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the term x 2 ^ in relation (JTJ. Near T c , at small x, 



m{x) = 1-2 e /3[(2VV^T7^)-i3.06] X 2 P (15) 

Here H = 0. That is, as the temperature approaches T c , 
m(:r) decreases with x more and more rapidly. 

Specific heat and susceptibility. — Substituting result 
(|14fl into formula for the free energy readily gives 
the specific heat, tC(T) = -Td 2 F/dT 2 . C(T > T c ) = 0, 
as is usual for mean- field theories. If T < T r , 



C(T) = 



(Tree) 2 



8(/3-l/4>' 



■ exp 



(16) 



where (nce) 2 /8 = 22.01 . . .. 

Similarly, one can consider the case of a non-zero ho- 
mogeneous magnetic field. Here we present the resulting 
expressions for the magnetic susceptibility: 

X (f3 > 1/4) = /r 1 - 1 , 



X (/3<l/4) = {l-y/T=4fl/(l + y/T=0).(l7) 

There is a finite jump of the susceptibility at the phase 
transition point: x\fi= (V 4 ) - ] =1 and x\P = (V 4 ) + ] = 3 
[see Fig. He)]. 

Response to a local magnetic field. — In networks, in- 
stead of correlations in space, one has to consider other 
options. In our case, to characterize the decrease of corre- 
lations between interacting spins, we use the distribution 
of linear responses to local magnetic fields. We apply a 
small field to the neighborhood of some point y: 

H(x, y) = h[9(x -(y- A/2)) - 6(x - (y + A/2))] (18) 

[9{x) is the theta-function, h and A are small] and find 
the change of the full magnetization which is induced by 
this field: fi(y) — J dx [m{x,y) — m(x)]. Knowing fi(y) 
readily gives the distribution P([i) of the response. 

Calculations are especially simple at T > T c . One can 
find m(x,y) by iterating Eq. @, which gives 



1 + VT - 4/3 



y 



-(l-VT=W)/2 



(19) 



Note the power-law divergence of the linear response as 
y — > 0. Relation l|19|) results in the power-law response 
distribution: 



PQi) oc ^-[1+2/(1-^1=33)] . 



(20) 



It is important that in contrast to "normal" continuous 
phase transitions, P{n) is a power-law function in the en- 
tire phase and not only at T c . As is natural, in the other 
phase, this distribution is a rapidly decreasing function. 
At the phase transition point, P(/i, (3 = 1/4) oc /i~ 3 . 



Discussion and conclusions. — Several points must be 
emphasised. 

(i) The spins on the oldest vertices are oriented in most 
of situations: m(x — 0) — 1 even above T Cl if any non- 
zero (positive) magnetic field is applied at least to one 
spin of the system. 

(ii) We considered networks with a specific annealing. 
The problem of a quenched disorder is more complex. 
However, there is a quenched situation, to which our re- 
sults are applicable. Let each new vertex have a large 
number N — greater than the final size of the network 
or of this order — of new connections to randomly chosen 
vertices. Let each of this edges bear the Ising coupling 
equal to 1/N. Then we arrive at the situation similar to 
that is shown in Fig. ^ 

(iii) The phase transition found in this paper, as well 
as the structural transition considered in Refs. 0, 0, 
13 El HIM El EHH, seriously differs from the usual 
BKT transition. In our case, the analogue of the power- 
law correlations takes place in the normal phase. In con- 
trast, in the traditional BKT transition, the power-law 
decay of correlations is in the phase with a non-zero order 
parameter. 

(iv) We stress that the more traditional-looking tran- 
sitions of Refs. [|| 0, H, IS LUJLL-3 are realized in equilib- 
rium networks where all vertices are statistically equiv- 
alent. In contrast, the networks, where we observed the 
transition with the BKT singularity, are specifically in- 
homogeneous. For our general conclusions, the specific 
1/j form of the inhomogeneity of the interaction in the 
Hamiltonian Q is important only in the region of rela- 
tively small j. Deviations from this form at larger j do 
not change the critical behavior. We studied a growing 
network but the problem has been reduced to the Ising 
model on a compact system with strong long-range in- 
homogeneity. We believe that our results are applicable 
to other systems with inhomogeneity of this kind. Fur- 
thermore, the Ising model is only a simple example of 
cooperative models were the observed transition should 
be present. 

In conclusion, we have solved the ferromagnetic Ising 
model on a highly inhomogeneous growing net. In this 
system we have found an infinite order phase transition 
with the BKT singularity. This transition separates a 
phase, where the distribution of the response to a local 
field is power-law, and the phase, where this distribution 
is rapidly decreasing. We suggest that this transition also 
occurs in other cooperative models on compact substrates 
with strong long-range inhomogeneity. 
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vice de Physique Theorique, CEA-Saclay for hospitality. 
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